Phase diagram of the extended Hubbard chain with charge-dipole interactions. 
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We consider a modified extended Hubbard model (EHM) which, in addition to the on-site 
repulsion U and nearest-neighbor repulsion V, includes polarization effects in second-order per- 
turbation theory. The model is equivalent to an EHM with renormalized U plus a next-nearest- 
neighbor repulsion term. Using a method based on topological quantum numbers (charge and 
spin Berry phases), we generalize to finite hopping t the quantum phase diagram in one dimen- 
sion constructed by van den Brink et al. (Phys. Rev. Lett. 75, 4658 (1995)). At hopping 
t = there are two charge density-wave phases, one spin density-wave phase and one interme- 
diate phase with charge and spin ordering, depending on the parameter values. At t the 
nature of each phase is confirmed by studying correlation functions. However, in addition to 
the strong-coupling phases, a small region with bond ordering appears. The region occupied 
by the intermediate phase first increases and then decreases with increasing t, until it finally 
disappears for t of the order but larger than U. For small t, the topological transitions agree 
with the results of second order perturbation theory. 
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I. INTRODUCTION 

Atomic polarizabilities are very important in solid 
state physics, since they contribute to the screening 
of the effective on-site repulsion. For example, from 
atomic data the bare on-site Coulomb energy for Cu"*"^ 
is Uo = E{Gu+^) + E{Cu+) - 2E{Cn+'^) ~ 17 eV, 
while constrained densitys- functional calculations for 
La2Cu04 give f/ - 10 e\WE. In the case of BaBiOs, 
while C/o(Bi+^) ~ 11 eV, it has Jpeen suggested that 
the effective U becomes negativeH. Interesting effects 
of the polarizability also appear in the physics of small 
clusters. In particular, while in neutral clusters of a few 
atoms of Hg the only relevant interatomic interactions 
are van der Waals forces (originated by fluctuations of 
the polarization) , a transition to covalent and finally to 
metallic bonding takes place as the size of the cluster 
increascsBij. A very good fit of the observed ioniza- 
tion potential as a function of cluster sizeQ has been ob- 
tained using a theory which includes charge-dipole and 
dipole-dipole interactions in a two band model (repre- 
senting the 6s and 5p electrons of Hg) .eI 

Recently, a simple one-band model includinjS .screen- 
ing effects has been proposed and studied.ljM The 
main simplifying assumption is that the electronic 
transitions responsible for the screening do not in- 
volve the valence electrons (only transitions from one 
deep level to another one well above the Fermi energy 
are considered). Furthermore, only nearest-neighbor 
charge-dipole interactions are included. These facts 
renders the model more amenable to many-body treat- 
ments, and interesting conclusions have been ex- 
tracted. In particular, in one dimension (ID), in addi- 
tion to the usual charge density wave (CDWl) and spin 
density wave (SDW) phases of the extended Hubbard 



chain, two neW||phases appear in the ground state for 
hopping t — a charge density wave with a unit 

cell of four lattice parameters (CDW2), and a com- 
bined spin and charge density wave with a unit cell 
of six lattice parameters called intermediate (I) phase. 
The ID phase diagram for i = was constructed, and 
the energy gap, static correlation functions and one- 
particle excitation spectra havftpbeen studied numeri- 
cally for a few cases with t ^ 0.BI13 However, the phase 
diagram for t ^ has not been investigated yet. 

The main difficulty in obtaining the phase diagram 
for t ^ is that, for finite systems, the conventional 
order parameters are smooth functions of the param- 
eters of the model, and the values at which the tran- 
sitions take place in the thermodynamic limit cannot 
be determined precisely. However, recently it has been 
shown that in certain cases, including the usual ex- 
tended Hubbard model with repulsive interactionsLj, 
the thermodynamic phases are characterized by dif- 
ferent topological numbers (determined by charge and 
spin Berry rakas^), and these jump suddenly for any 
system size.lljil3 Extrapolating the position of this 
jump to the thermodynamic limit, a very accurate de- 
termination of phase boundaries was obtained, as es- 
tablished by comparison with available exact results 
in the strong-coupling lirpit and field-theory ones in 
the weak-coupling caseOEj In ID, the method of the 
Berry phases is related to the method of crossings of 
excitacion levels based on the confojarpal field theory 
with renormalization group analysis.E§E3 This method 
(abbreviated as level-crossing method in the rest of the 
paper) is not expected to work in the strong coupling 
limit. 

In this work we use the method of Berry phases to 
calculate the phase diagram of the model for t ^ 0. 



In particular, we show that in this case all the ther- 
modynamic phases can be characterized by the topo- 
logical numbers, including also a new phase with_bond 
ordering previously suggested in the literature.E£l The 
nature of each phase is confirmed by calculating differ- 
ent correlation functions. In Section II we describe the 
model and the strong-coupling thermodynamic phases. 
Section III contains a brief description of the calcula- 
tion of the topological numbers and the values they 
take for the t = phases. In Section IV we present 
the quantum phase diagram for several values of Uq /i, 
compare it with the results of perturbation theory in t, 
and discuss the behavior of the correlation functions. 
In Section V we summarize our results and discuss the 
extent and limitations of the method of the topological 
transitions and that of the level crossings. 



II. MODEL AND T = GROUND-STATE 
PHASES 

The model is defined by the Hamiltonian 

H = Ho + H,d, (1) 
where is the usual extended Hubbard chain: 

H =-t ^(cj+i^c,^ + H.c.) + C/o XI "''T^a 

i 

(2) 



and Hcd describes charge-dipole interactions (dipole- 
dipole interactions are neglected): 



Hcd — — X Pi ■ Ei. 



(3) 



Here Pi and are the dipole moment operator and 
electric field at site i respectively. Treating this term 
in second-order perturbation theoryOild it takes the 
form — i aiEf , where ai is the polarizability at site 
i. Considering only the contribution to the electric 
field coming from nearest-neighbor atoms, and assum- 
ing that ai is independent of the site and its occupancy, 
one can write: 



Hcd = ~P'^{ni+i - n^-i) 



(4) 



where P — ae'^/{2a'^) and a is the lattice parameter. 
Expanding Hcd, the total Hamiltonian can be cast into 
the form: 

^ = ^tJ2i4+l.c^. + H.C.) + (Uo - 4P) J2 

ia i 

+ vY,n^n,+l + 2PY,n^n^+2~2PJ2n^. (5) 



The last term merely rftnarmalizes the chemical po- 
tential. As stated beforeErll3, within the above men- 
tioned approximations Hcd reduces the on-site interac- 
tion and generates a next-nearest-neighbor repulsion. 



We are interested in the half-filled case of one par- 
ticle peij-site. For t — the model can be solved 
exactlyEni^, and at this filling four different phases ap- 
pear in the ground state depending on the parameter 
values. We use the notation Q, and |, to rep- 
resent empty, doubly occupied, and singly occupied 
sites with spin up and down respectively. If V is much 
larger than the other interactions, nearest-neighbor oc- 
cupancy must be avoided and the ground state can be 
represented as (... O • O • O • O •■■■)• This state is 
denoted as charge density wave 1 (CDWl). If Uq is the 
dominant interaction all sites are singly occupied, and 
an infinitesimally small t is enough to introduce anti- 
ferromagnetic correlations between nearest-neighbors. 
Thus, we represent this spin density wave (SDW) as 
(•■• Tititi ••■)i although there is no true spin long- 
range order in one dimension. If P is very large, oc- 
cupancy at next nearest neighbors is avoided and the 
ground state (CDW2) is (...QO ••00 •••■•)■ Finally, 
for intermediate values of the parameters (for example 
V — O.bUo, P = O.lbUo) there is an intermediate (I) 
phase, the unit cell of which can be represented as • t 

O^iO- 



III. BERRY PHASES 

We have calculated the ground state Igi^i, $x)} 
of the model in rings of different lengths L, threaded 
by fluxes $cr for spin a (the hopping for each spin 
is affected-|by, a factor e**"/^ as described in detail 
elsewhereliirEJ) . We call charge (spin) Berry phase 7c 
(7s), the phase captured by the ground state when it 
is followed adiabatically in the cycle < $ < 2tt, keep- 



ing $1 = $1 = $ (<I)| 



-$l = $). Discretizing the 



interval < $ < 27r into N points $i = 2'Ki/N, the 
Berry phases are calculated numerically by the gauge- 
invariant expressiorU: 



7c 



hm {Im[ln(H^-2(5(<i>„±cO,)|g(<I,,+i,±<i>,+i)) 

A/ — i-oo 



(g($w_i,±$^_i)|5(27r)))]}. 



(6) 



where \g{2TT)) represents \g{2TT,± 2tt)) obtained di- 
rectly from \g{0, 0)): 

|5(27r)) ^expiz^^j-e^:,! ±'^.i)]l.9(0,0)). (7) 



An important property of the charge Berry phase is 
that if the system is modified by some perturbation, 
the change in the polarization P| is proportional 
to the corresponding change in 7c.c-j Here Pg. is the 
contribution of electrons with spin a to the polariza- 
tion of the system. Similarly, changes in 7^ are related 
to changes in the difference P^ — P| between the electric 
polarizabilities for spin up and downE^: 



AP| ± APi = eA7c,s/27r [mod(e)]. 



(8) 



A crucial property for our purposes is that, in sys- 
tems with inversion symmetry like ours, 7c and 7^ can 
only be either or tt [mod(27r)] (the argument of the 



logarithm in Eq. m) turns into its complex conju- 
gate under inversion). Then, the Berry phase vector 
7 = (7c, 7s) cannot vary continuously, and a jump in 
7 corresponds to a transition in at least one of the 
topological numbers 7c/7r, 7s/7r. Our main hypothesis 
is that the topological transitions correspond to phase 
transitions in the thermodynamic limit. This has been 
confirmed by different numerical and analytical .st.i i d ies 
in the Hubbard model with correlated hopping.EjO In 
the present case, we check this by comparing different 
correlation functions (see next section). In addition, 
it has been shown that the opening of a spin gap in 
a Luttinger liquid phase of a spin rotational invariant 
model is accompanied by a topological transition in 

Let us discuss the values of 7 for the four possible 
ground states of Eq.(|^) in the half-filled system with 
t = 0, discussed at the end of the previous section. 
For t = 0, we can choose a gauge in which all scalar 
products in Eq. (^) are equal to 1, except eventually 
the last one. Then the Berry phases are defined by 
the argument of the exponential in Eq. (^. If the 
ground state is the state (... O • O • O • O ••■•) of the 
CDWl, then this argument is clearly for 7^, while for 



7c it is equal to i'^ j("iT +'^ji) = «T Em=i 4"i = 
i2iT{L+2) = [mod(27rj)]. The same result is obtained 
for the other CDWl state [... • Q • O • O • 0-], 
so 7(CDW1)= (0,0). One can construct a Neel state 
(representing the SDW) displacing half of the electrons 
(for example those with spin up) one lattice parameter 
either to the left or to the right. Then, from Eqs. (^ 
one has 7(SDW)= (tt, it). This result can also be easily 
obtained calculating the expopant of (|^), and is not 
altered by spin fluctuationsO'EJ For the CDW2 with 
order (... O O • • 00 • • •■•) clearly 7^ = 0, because 
the spin at each site is compensated. For 7c one has 
xE,j(n,T +".i) = xEtC\(16m-2) = 2^(L + 
4- 1/2) = TT [mod(27r)]. Then, 7(CDW2)= (7r,0). 
While the phases above have been calculated for t = 
0, it is clear that by continuity each thermodynamic 
phase retain the same value of 7 unless a topological 
transition takes place. 

The vector 7 in the I phase cannot be calculated 
as easily as above for i = 0. The reason is that the 
state (...• t O* I O---) does not have a definite parity 
under inversion and other states related with it by in- 
version have different exponents of Eq. (|^) . For t ^ 
the ground state should have a well defined parity, and 
7 depends on the particular linear combination of the 
above mentioned states which constitutes the ground 
state. Numerically we found 7(1)= (0,7r). This is a 
fortunate fact, since then at t = the four thermo- 
dynamic phases are characterized by all the possible 
values of 7. Thus, the corresponding phase bound- 
aries can be determined from topological transitions. 
Notice, however, that this does not rule out the pos- 
sibility that at i 7^ new phases having the same 7 
values could appear. We will show in the next sec- 
tion that this is in fact what happens. To summarize 
and for future use, the Berry phase vectors for each 
of the strong-coupling phases are 7(CDW1)= (0,0), 



7(SDW)= (7r,7r), 7(CDW2)= (7r,0), and 7(1)= (0,7r). 



IV. NUMERICAL RESULTS 

In this section we present the ground-state phase 
diagram of model (||) for t ^ 0. As discussed in the 
previous section, it has been obtained by studying the 
behavior of the Berry phases in rings of L = 8 and 
12 sites. Notice however that the order corresponding 
to the I phase cannot be accommodated in rings of 8 
sites because of its 6-site unit cell. In Fig. 1 we show 
the phase diagram for Uo/t = 10, 3 and 1, which cor- 
respond to the strong to moderate coupling regimes. 
As can be seen in Fig. la, for Uo/t — 10 the results 
are similar to the strong limit case {t — 0). Moreover, 
there are no sensible differences between the values ob- 
tained in rings of 8 and 12 sites, which indicates that 
they are representative of the thermodynamic limit re- 
sults. Fig. lb shows that already for Uo/t — 3 the 
behavior departs from the one corresponding to the 
strong coupling limit; in particular, there are stronger 
finite-size effects in the transition lines CDWl-I and 
CDW2-I, especially when the transition line is deter- 
mined through the jump in the spin Berry phase. The 
line CDW2-I is better determined by the changes in 
the charge Berry phase, which shows no major finite- 
size effects. We expect that, for this line, going to 
larger systems the jumps in the spin Berry phase will 
coincide with those of the charge Berry phase. 

Note that with a small increase in t from the strong 
coupling limit, the I phase first increases with respect 
to the CDWl and CDW2 phases. This can be under- 
stood from perturbation theory in t. The energies of 
the different phases calculated up to second order in t 
are the following: 



4t^ In 2 



E{SDW) = N [V + 2P 



Uo-^P-V 



' J J 9/ 
i.(CDWl) = A. + 2P-^^ 



i^(CDW2).A. ['^-^p+y- ,,p!y_u^ 



Uo + 2V -t^ [- 



1 



P 8P + V -Uo 



(9) 



The phase diagram that results from these energies is 
represented in Fig. 2. For Uo/t = 3, the second or- 
der expression for £'(SDW) is no more valid (it gives a 
very low energy). Thus, we have not drawn the bound- 
aries of the SDW phase for this value oiUo/t. Instead, 
the other boundaries of the I phase, and in general all 
boundaries for Uo/t = 10 (except when some denomi- 
nator becomes of the order of t), agree very well with 
the topological transitions represented in Fig. 1. 
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FIG. 1. Quantum phase diagram of model (|l|). a) 
Volt = 10; b) Uo/t = 3; c) Uo/t = 1. Solid dots and 
triangles correspond to jumps in the charge Berry phase 
obtained by considering rings with L = 8 and 12 sites re- 
spectively; open dots and triangles indicate jumps in the 
spin Berry phase for systems of the same length. The thin 
lines correspond to the t = phase diagram. 

From Eqs. and Fig. 2, we see that near the 
points where the CDWl, CDW2 and I phases meet for 
t = the stabihty of the I phase increases with t. This 
is due to the fact that the electrons of the doubly oc- 
cupied sites of the I phase can jump to the neighboring 
singly occupied sites with small energy cost in P, and 
without increasing on-site and nearest-neighbor repul- 
sion energies. As a consequence, for small t ^ 0, the 
electrons in the I phase become more mobile and gain 
more kinetic energy than in the other two phases. 

On the other hand, for Uo/t — 3 the topological 
transitions seem to indicate that a new phase appears 
between the SDW and CDWl phases (see Fig. 1). 
This was only slightly visible for Uo/t — 10 and is out 



of the reach of the perturbative calculations. We wave 
studied this possibility in more detail by considering 
the scaling of critical points on the axis P = shown 
in Fig. 3. From this figure we see that the jumps in the 
spin and charge Berry phases will not coincide even in 
the thermodynamic limit; this fact is pointing to the 
emergency of a new order termed bond ordering wave 
(BOW), as suggestecU|n a recent study of the phase 
diagram of the EHM.E3 
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FIG. 2. Phase boundaries obtained by perturbation 
theory in t. Full line f/o/t = 10, dashed li«e,;7o/f = 3. The 
dotted line is the known result for t = oJatJ. 

For U{)lt — 1 (see Fig. Ic) the phase diagram 
looks quite different to the strong coupling case, with 
the I phase shrinking from a triangular-like shape to 
a smaller biconvex region. Furthermore, the BOW 
phase remains present without much change. Notice 
that using bosonization (valid in the extreme weak- 
coupling xegime i ^ oo), the BOW phase is absent for 
P = oJlaM and the I and CDW2 phases should also 
disappear since the expected instabilities have wave 
vector 2A;f = ti"- 
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FIG. 3. Values of Vjt at which the charge (solid dots) 
and spin (open dots) Berry phases jump on the P = axis, 
as a function of the ring length L. We have taken (Jo /t = 3. 

In order to have a better characterization of the or- 
der present in the different phases, we have also com- 
puted static correlation functions along the P — Q 
axis from Vjt = \ lo Vjt ~ 2, that is, travers- 
ing across the SDW, BOW and CDWl phases. 
Fig. 4 shows the spin and charge correlation func- 



tions at the largest distance L/2 = 5 for a 10- 
site chain, and also the correlation function corre- 
sponding to the expected order in the BOW phase, 

CBow(i/2) = ([Ob{L/2) ~Ob\[Ob{0) -Ob\), with 



Ob(j) = (-PE.(4 



H.c.)- The vertical dot- 



ted lines correspond to the transition points extrapo- 
lated in Fig. 2, and we can see the expected behavior 
for these three functions along the chosen path. 
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FIG. 4. Static charge (Cc), BOW (Cb) and spin 
(Cs correlation functions at the largest distance L/2 = 5 
for a 10-site ring. The dotted vertical lines indicate the 
extrapolated points in Fig. 3. 

Furthermore, in Fig. 5 we plot the fluctuations 
of the three order parameters [for example xbow = 
e'9(-R--R^)(OB(«)OB(i))] at their peak values 
(corresponding to fc = tt), as a function of the ring 
length and for V/t = 1.5, that is, inside the BOW 
phase. Both the spin and charge fluctuations dimin- 
ish with increasing lengths, while the BOW fluctua- 
tions increase and most likely dominate in the ther- 
modynamic limit. The same analysis in the SDW and 
CDWl phases is not so conclusive because the BOW 
fluctuations do not clearly decay when increasing the 
size of the system; however, since the orders present 
for small and large values of V/t are not a matter of 
discussion, this must be only due to the small size of 
the systems considered. 
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FIG. 5. Charge (open dots), spin (open squares) and 
BOW (solid dots) order-parameter fluctuations at fc = vr 
as a function of the ring length L. We have considered the 
parameters values P = 0, V/t = 1.5 and Uo/t = 3, which 
correspond to a point well inside the BOW phase. 



V. DISCUSSION 

We have studied the effect of electronic polarization 
P on the phase diagram of Xha ID extended Hub- 
bard model. Following Refs.Q'EEl, we have neglected 
the screening effect of valence electrons, althougtuiit 
was important to explain experiments in Hg clusteraZla. 
At P = 0, for a bipartite lattice in any dimension two 
phases are expected: a SDW for large on-site Coulomb 
repulsion Uq (without true long-range order in ID), 
and a CDWl for larg&*earest-neighbor repulsion V. In 
ID the results of Ref Jl3 and ours support the existence 
of a narrow region of long-range bond ordering between 
the other two. phases, if the hopping t ^ 0. Previous 
studies in IIJJ'E] have shown that for small t and suffi- 
ciently large P, two new phases appear: CDW2 and I, 
with unit cells containing four and six lattice parame- 
ters respectively, as described in section IL 

For i = it is clear (see section II) that a transition 
from one phase to the other is accompanied by charge 
and /or spin transport. Since the charge Berry ^flhase 
7c is a measure of the polarization of the system,E3 and 
similarly a jump in the spin Berry phase 7^ indicates a 
jump ia the difference in polarization for spin up and 
down,t3 it seems appropriate to determine the phase 
diagram froin the jumps in the topological numbers 
7c/7r, 7s/7r.ll3 Using this method, we find that for fixed 
Uq, V and P, the regions occupied by the CDW2 and 
particularly the I phase, are reduced for sufficiently 
large t. This result is expected on general grounds, 
since only instabilities with wave vector tt are expected 
in the weak coupling limit {t — > 00), which should .be 
well described by continuum limit field theory.oliSlij 
For small t (up to Uo/3) the region occupied by the 
I phase first increases with t, and this effect can be 
understood using second-order perturbation theory in 
t. 

Although in the present case we succeded in deter- 
mining all phase boundaries from the jumps in the 



topological numbers 7c/7r, 7s/7r, there are some cases 
in which the thermodynamic phases are character- 
ized by the same topological numbers, and the phase 
transition cannot be detected. In other cases some 
topological transitions do not have an obvious phys- 
ical meaning. t3'E3 An example of the former is the 
transition between SDW and triplet superconductiv- 
ity which takes place for positive U and negative V 
in the extended Hubbard model in ID, which is how-i 
ever well determined by the leiieLcrossing method.ll3 
This method, as used in Refs.ll3'E£l, is limited to ID 
and is based on the fact that for critical systems, due 
to conformal invariance, the exponent of the distance 
dependence of different correlation functions is given 
by the corresponding excitation energies £j For sys- 
tems at (or slightly out of) criticality, described by 
the sine- Gordon model at low energies, renormaliza- 
tion arguments were used and logarithmic finite-size 
corrections can be eliminatedaisiag-pppropriate combi- 
nation of excitation energies. ESEifEj This method can- 
not be applied to phases which are out of the reach of 
the continuum-limit field theory like the CDW2 and I 
phases, or when there is phase separation. Instead, in 
the strong-coupling limit {t 0) it is usually very sim- 
ple to calculate tlie_Berry phases, even in some phase 
separated statesE^o (it reduces to a simple counting 
argument, as shown in Section II), and the method of 
the topological numbers is very precise. It is inter- 
esting to note that, in spite of their different meaning 
and expected limitations, both methods give the same 
results for transitions among the SDW, CDWl and 
BOW phases. The method oLJ3erry phases has also 
been applied to 2D transitionsEZl. 
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